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ABSTRACT 


The increasing weight and speed of carrier based aircraft are 
ee the limit of conventional piston type energy absorbers which 
are used for arresting gear. Turbo-type absorbers have been proposed 
as an alternative, This study investigates a turbo-type energy 
absorber for an aircraft arresting gear under development by the 
Naval Air Engineering Center. A mathematical analysis, computer 
Pamilation and performance prediction are given for each mode of 
operation, It was initially expected that the absorber could operate 
in just two on forward and reverse, but the analysis shows that 
-four distinct modes are possible, since not only shaft rotation but 
the path of fluid Flom can be reversed. Theoretical performance 
predictions are also compared with test data from an existing smaller 
scale version of the sneeabes, The agreement is excellent. An 
approximation of scaling effects on power absorption is also 
presented. It is concluded that the turbo-type absorber is basically 
adequate to meet the demands of carrier operation for the forseeable 


future. 
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I. INTRODUCTION 


Over the years, the speed, size, and weight of carrier based air-— 
craft have steadily increased, thus demanding an ever larger energy 
absorption capacity from the shipboard arresting gear. Up to now, 
design improvements to meet these increasing demands have been limited 
to mere enlargement of the basic linear hydraulic energy absorber, 
which has been of the piston type. Modern aircraft, however, have 
loaded this system to its upper limit, and make any further development 
along these lines impractical due to weight and space limitations. The 
need exists for a new concept which will minimize space and weight 
demands while providing the necessary increase in energy absorption 
capacity. 

The solution proposed involves the use of a turbo-type rotary 
hydraulic energy absorber. In essence the fluid is put into motion by 
a suitably designed rotor and is wedda by a stator and by the 
resistance of the flow passage. The power input to the rotor is 
absorbed by the hydraulic resistance of the system, and is eventually 
dissipated as heat. 

The energy absorber is directly coupled to a cable storage drum. 
The drum profile is programmed to match drum rotational velocity to that 
required for optimum retarder performance. A device of this type was 
introduced in 1968, and is known as the M-21 energy absorber. The 
initial design gives satisfactory performance, but is of rather limited 
capacity, thirty million ft.lbs. Therefore, a program has been under- 
taken to redesign the absorber to increase the capacity to fifty-five 


mation ft. lbs. 
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A development program was initiated by the Naval Air Engineering 
Center (NAEC) (Ref. 1). A full scale test of the programmed drum 
coupled with the M-21 absorber was conducted to determine the 
feasibility of the overall concept and to obtain performance data 
needed in the redesign program. A schematic diagram of the proposed 
system is shown in Figure 1. 

A theoretical approximation of the flow conditions in a machine 
of this type was derived by Professor Michael H. Vavra (Ref. 2a te 
facilitate investigation of the effects of various design parameters 
NAEC formulated a computer program. From the computer analysis, a 
workable er eee of design parameters was chosen such that torque 
and energy requirements were satisfied at the design point, the con-— 
dition at which the absorber is expected to perform most often. 

Professor Vavra has also been commissioned to develop and test a 
flow model of the new design. This model will be useful for refining 
losses, checking for flow separation, and obtaining information that 
is not available from the mean streamline analysis. 

In this paper, a theoretical simulation and performance analysis 
of the energy absorber in all modes of operation is presented in terms 
of dimensionless parameters that facilitate comparisons of various 
designs and interpretation of model tests, The retraction problem 
is also studied, and preliminary estimates of the power needed to 
retract the programmed drum are given. 

A throttling system is incorporated into the design to aid in 
arresting aircraft of various sizes. The discussion and analysis 


presented here are concerned with the condition of open throttle only. 
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II. DESIGN VALUES AND CONSIDERATIONS 


The main features of the energy absorber design that is analyzed 
here are summarized in Figures 2 and 3, and Table I. 

The customary approach in the design of turbomachinery is first to 
optimize conditions with respect to a single specified design point. 
The resulting performance over a broad range of off-design conditions 
is then investigated to determine overall acceptability. The condition 
at which the absorber is expected to operate most often is naturally 
chosen as the design point. Specifically, a design torque of 
340,000 ft.1bs. eewaed at a rotational speed of 76 rad./sec. 

-The maximum rotational speed the cable will impart to the drun, and 
consequently to the rotor, is about 83 rad./sec. This limit is fixed 
by the allowable landing speed of ime aircraft, and by drum design. 

The mechanical power imparted to the rotor is continuously 
converted to thermal energy of the fluid. The temperature rise of the 
fluid is significant and imposes a constraint on the design. The 
temperature limit is the maximum allowable temperature the fluid may 
attain without causing cavitation. Obviously, heat must be removed to 
Maintain equilibrium conditions. For this reason, provisions are made 
to bleed off a portion of the circulating fluid continuously and to 
replace it by fresh cool fluid. During a rapid series of arrests, the 
cooling problem is especially critical. 

The fluid used in the absorber is an ethyl glycol mixture. This 
fluid is acceptable because it remains usable over a wide range of 


temperatures. 
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Figure 3 DIAGRAM OF ABSORBER BLADING 


Ay 


Table I PARAMETERS OF THE AS3SORBER 


R, = 13.61 inches 
R, = 16.25 inches 
bs =) 16.07 aianches 
A =l14,.930inehes 
bp = 2.02 inches 
ba = 2.02 inches 
are= 4.18 inches 
r= 2.02 inches 
(L a iy 
= 10.26 inches 


(hy + ty) 


‘Number of blades for rotor and stators [72 
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Fluid properties: 


P2095 Slues/fte cp = 016 


On) 
i] 


i] 
. 
—y 


Fai fo 407? Ft. “/sec c 


< 
N 


18 


A detailed discussion of the effects of various parameters and loss 
coefficients is undertaken in the analysis section. A few words, 
however, are in order here. The radius of curvature rs of the inner 
contour of both the upper and lower passage must not be smaller than 
the blade height. This insures that the flow does not separate, and 
it also increases absorption capacity. Estimates of losses in this 
machine are obtained through various empirical turbomachinery cor- 
relations as is shown in the analysis. Enthalpy losses in the passages 
are significant. Much of this loss is due to the bending of the 
passage with friction responsible for the remainder. However, a 


change in the coefficient of friction c, which expresses the direct 


1g 
wall friction losses in the passages has little influence on overall 


performance. A value of c, = 0.016 was chosen since it gives best 


f 
agreement with M-21 test data. 

The blade shapes first optimized by Professor Vavra from a 
performance standpoint were found by NAEC to be unacceptable ona 
stress basis. Blade thickness, spacing and angles were accordingly 
revised by NAEC, and a blade of sufficient strength was found; its 
characteristics are listed in Table l. | | 

The reliability of this analysis can only be judged experimentally 
by correlation with test data. The theory correlates very well with 
actual test results on the M-21 design. From this it appears that the 
theoretical analysis is reasonable and that the new design should be 
workable, if not optimum. 

fered Ga operations often require a series of rapid and 


repeated arrests. The drum must be retracted thus reversing the 
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rotational direction of the rotor. This process is investigated to 
obtain preliminary estimates of the power required to retract the 


absorber. 
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III. ANALYSIS OF THE PROBLEM 


A. MODES OF OPERATION 
The energy absorber is capable of operation in four modes, These 


modes are summarized in Table II. 


Table II MODES OF OPERATION 


Flow Direction Shaft Rotation 
Normal Reverse 
Normal Mode 1 Mode 2 
Reverse Mode 3 Mode 4 


Mode 1 corresponds to design point conditions. In this mode, the 


A A 
L} 


fluid passes successively tnrough stations 1, 2, 3, 4 and back 
in that order. Note that the flow leaving the rotor at station 2 is 
moving radially outward. See Figures 2 and 3. For fluid circulation 
in this sense, conditions at stations 2 and 4 are particularly sig- 
nificant. At the rotor exit, station 2, the direction of the flow 
with respect to the rotor, is assumed to remain constant at all flow 
rates. At stator exit, station 4, the direction of the flow with 
respect to the stator is likewise dssumed to sent constant. While 
this assumption of fixed relative and absolute flow directions is not 
exact, it represents an excellent approximation when the blades are 
closely spaced as in the present design. 


In Mode 2, the sense of the fluid circulation is exactly as 


described above. However, the direction of shaft rotation is reversed. 
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Naturally, stations 2 and 4 are still governing in the sense that flow 
directions at these stations are known, 

In Mode 3, the shaft rotation is the same as for Mode 1, but the 
direction of fluid circulation is reversed, Consequently, the fluid 
passes successively through stations 2, 1, 4, 3 and back to 1, in 
that order, Note that the rotor exit is now at station 1, and the 
direction of flow at this station is radially inward. In this case, 
stations 1 and 3, are controlling. The relative flow direction at 
rotor exit, station 1, and the absolute flow direction at stator exit, 
station 3, are now assumed to remain fixed. 

In Mode 4, the sense of fluid circulation is the same as in 
Mode 3, but the direction of shaft rotation is reversed. Of course, 
stations 1 and 3 still govern flow direction as in the previous case. 

The same set of equations describe Modes 1 and 2, but a separate 
set of equations is needed to describe the reverse flow Modes, 3 and 4. 
Note that only two sets of equations are needed for the four Modes, 
and each set is governed by the flow direction relative to tk blade 
exits. Two separate computer programs are needed since one program 
will be used to solve each set of equations. Program Normal finds the 
solution for Modes 1 and 2, and ieee Reverse finds the solution for 
Modes 3 and 4, Both programs will be discussed in the appendix. All 
quantities and results discussed throughout the paper are obtained 
from these programs. 

The velocity diagrams for each mode are presented in Figures 4, 5, 
a and es All angles and velocities are drawn to scale. 

The thermodynamic processes in the absorber are shown on the 


enthalpy-entropy diagrams, Figures 8 and of The process for normal 
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Figure 4 Velocity Diagrams for Mode 1 Operation at 
2 

Ww. R Z 
= ue Gees = 4.8 x 40° 


Ne 


v 


a3 


U2 
2 
W 
ze 
ee 
U; 


ep ea 


Scale; O.5 inch= 25ftAec 


Figure 5 Velocity Diagrams for Mode 2 Operation at 
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scale: 0.5 inch=100 ft/sec 


Figure 6 Velocity Diagrams for Mode 4 Operation at 
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Figure 7 Velocity Diagrams for Mode 4 Operation at 
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circulation, Modes 1 or 2, is shown in Figure 8, while that for 
reverse circulation, Modes 3 or 4, is shown in Figure 9. The 
direction of shaft rotation affects the exact quantitative enthalpy 
values. The graphs, however, are merely qualitative. 

A "total head" H (per unit mass) is used which differs from the 
standard total enthalpy Hne This is discussed in greater depth in the 
analysis section, but a few words are needed at this point to clarify 
the discussion, 

The fluid is considered incompressible which allows the thermal 
and mechanical effects to be uncoupled. A "total head" is 


conveniently defined as 


P 2 
race ee) ag 
p p 2 


Henceforth H is referred to simply as head. The head and total 


‘enthalpy are linked by the simple relation 


=f > 2 


where u represents internal energy per unit mass. 


B, DIMENSIONAL ANALYSIS 
Results are presented in terms of appropriate dimensionless 
quantities. The use of non-dimensional terms simplifies the cor- 
relation of data, the comparison of various designs, and the inter- 
pretation of model tests. 
_ The three quantities chosen as consistent dimensional reference 
parameters are the radius of the rotor outlet Ros the fluid density Py 


and the shaft rotational speed w or Ovrax? whichever is appropriate. 
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Using these parameters, one defines the effective Reynolds number 


for this design to be 


where v is the kinematic viscosity. And the dimensionless total head is 


In the discussion of results, it is convenient to define the flow 


coefficient ¢ in the following specific way, namely, 


sg a2 
e wR, 


where the algebraic signs of w and $ are always positive for normal 
shaft rotation and negative for reversed shaft rotation. This sign 
convention for shaft rotation is maintained regardless of the direction 
of the fluid circulation, hence the use of the absolute magnitude es 
in the definition. 

In the actual computer programs NORMAL and REVERSE, however, some 
what different sign conventions are more convenient. All equations are 
written in such a way that the three quantities Vigo WR, and ¢ are always 
positive. For this reason the four modes are computed as separate cases. 

For comparison of various modes and designs the dimensionless torque, 
power and rotational Had coefficients are respectively 
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The detailed analysis considered later also useg§ various 
dimensionless quantities for convenience in the derivations. The 
parameters used to facilitate the non-dimensionalization in that case 
are not necessarily the same as introduced in the present section, and 


the separate cases should not be confused, 


C. ENERGY BALANCE. AT EQUILIBRIUM 
In general the rotor adds energy to the fluid. In traveling 
“through the passages and returning to its initial state the fluid 
dissipates this mechanical energy through various losses. Under 
equilibrium these losses exactly cancel the original energy input. 
It turns out that the energy input of the rotor can be expressed 


in terms of the flow coefficient $ in the form 


= $ 
4 Bay A+B 


where A and B are true constants, The losses can be expressed as 
aN a 4g Bb acc’ 4° 
“Hoss 7 


where al Bey and Gc may be considered constant to a first approximation, 
Actually these quantities are weak functions of 6, a fact which must 

be considered in the exact solution. ‘The energy relations whieh £ix 

the equilibrium state are shown in Figures 10 and 11. Figure 10 


represents the two modes which correspond to normal flow direction. 
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Figure 11 shows the two modes for reverse flow. These results are 
drawn to scale, except that the parabolic lines are only approximate. 
The final equilibrium points, where the losses balance the energy 
input, however, and the relative distribution of losses at these 
points are in correct scale. 

Initially, it was thought that only two of the above four modes 
would prove to be physically possible, one corresponding to normal 
shaft rotation, the other corresponding to retraction. Surprisingly 
enough, the present analysis has demonstrated that operation in all 
four of the modes is theoretically possible, This point is clearly 


indicated in Figures 10 and 11. 


D., MODEL TESTING 

Professor Vavra has proposed that a model testing program be under= 
taken. The model will not be a true dynamic model, but will be a purely 
static device that uses air as the working fluid. Construction of the 
model will constrain the air flow to the normal direction only, However, 
the scale model should be useful in checking for flow separation 
through the passages and refining the loss coefficients for the 
analysis of the absorber. The effects of the tangential component of 
velocity will also be investigated. 

The discussion of the above static model is not within the scope 
of this thesis. However, a hypothetical model is considered which is 
a true dynamic model of the actual full scale device. The size of the 
hypothetical model is taken as identical to that of the static model. 
The model is scaled to 0.569 of the size of the actual absorber. In 


both cases the fluid is air. When expressed in dimensionless terms, 
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the performance of the hypothetical model is basically the same as that 
of the full scale machine, except for Reynolds number effects. The 
full scale machine operates at a maximum Reynolds number, as defined 
above Of Al86. = 10° while the maximum Reynolds number of the model is 
56 (0.x 10°. The results for this model are presented along with the 


results for the real absorber, 


E, ASSUMPTIONS AND LIMITATIONS OF THE ANALYSIS 

A one-dimensional mean streamline analysis is obviously only a 
crude approximation of the actual performance of any turbomachine, 
but is a starting point necessary in most development programs. The 
validity of the assumptions used are very important to the reliability 

of the analysis. The primary assumptions used in this report are 
Pee ied below. 

The fluid is treated as incompressible, This assumption is 
excellent as long as the provisions to maintian the fluid at a 
temperature lower than ae cavitation temperature are adequate. This 
conveniently uncouples mechanical and thermal effects and allows a 
form of "total head" H to be defined which is independent of 
internal energy. 

If the velocity distribution across every section were turly 
uniform, the meridional component V. at. the mean streamline could be 


H 
computed simply from a continuity relation of the form 


where Q is the volumetric flow rate and Ay is the meridional cross- 


section area. The actual velocity distribution, however, is non- 


DD 
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uniform because of the effect of streamline curvature, boundary layers 
and so on. To allow for these effects, approximately, the actual 
meridional velocity Vy is assumed to be governed by a relation of the 


forn 
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where. the so-called blockage factor K, is treated as a constant. The 
blockage factors were extrapolated from M-21 test data. Their reliability 
in the reverse flow modes is doubtful, however, since all test data 

were obtained only from Mode 1 operation. An extension of these data 

to all modes would be very useful. 

The present analysis treats operation of the absorber as a steady 
state phenomenon, Strictly speaking, the operation is actually 
unsteady. However, the arrest and retraction processes require a 
paiiiiie tent number of drum revolutions that the steady state assumption 
provides a fair approximation. 

A major assumption is that the relative flow direction at the rotor 
outlet and the absolute flow direction at the stator outlet are 
constant. For the case of normal flow this is an excellent ap- 
proximation, However, the validity is reduced for the reverse flow 
modes. In Modes 3 and 4, the flow encounters a sharp leading edge, 
and exits past a blunt trailing edge. The implication is that the 
flow is not as well guided as in the normal flow case. Thus the 
actual turning angle of the flow is undoubtedly less than predicted 
by the theoretical analysis. The power absorbed in both reverse flow 


modes would consequently be less than predicted. 
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Losses, in general, are approximated by various empirical cor- 
relations and are reasonable as shown by test results. They are 
discussed in depth in the detailed analysis section. 

Incidence losses, however, pose a problem. The incidence loss 
approximations introduced in the analysis section are reasonable as 
long as the incidence angles remain small. This is obviously not 
the case for the retraction modes as shown in Figures 5 and 7. The 
approximations break down physically at these severe incidence angles. 

It was thought earlier that two of the four theoretical modes 
might show themselves to be physically impossible. For example, a 
large incidence loss might conceivably preclude the attainment of the 
necessary energy balance. This has not proved to be the case, however, 
It appears from the present analysis that operation in all four modes 
remains possible physically. It is true, of course, that the predicted 
performance for the two retraction modes, Modes 2 and 4, will be 
inaccurate because of these large incidence effects, but this in itself 
does not eliminate the existance of these modes. 

At extreme angles of incidence, a distinct possibility of stalling 
exists. In normal operation, this limits the attainable power absorption, 
In a retraction mode, however, stalling would not be totally undesirable. 
If the rotor stalls, the energy absorbed by the device is greatly 
reduced, and the power necessary to retract the rotor is minimal, 

All of the foregoing assumtions are thought to be reasonable. 
Refinements of various points would be useful, especially for large 


incidence angles and blunt trailing edges. 


IV. RESULTS AND DISCUSSION 


The discussion of the problem up to this point has already included 
a Piet of significant results reached in the analysis as shown in 
Figures 4 through 11. 

Additional detailed results for all modes of the absorber are 
presented in Figures 12, 15, 14, and in Table III. The results for 
Mode 1 of the dynamic model are also included. 

The main results of this sree tien are summarized concisely in 
Figure 12. This is a logarithmic plot of the dimensionless power II 


as a function of the dimensionless shaft speed w/w, Performance 


MAX* 
-of the proposed absorber, of the hypothetical model, and of the 
existing M-21 machine are compared on this graph. Note that the 
theoretically predicted performance of the M-—21 version as shown by 
the continuous line is a excellent agreement with the actual measured 
performance as shown by the data points. 

The actual horsepower is plotted as a function of absolute shaft 
speed for all four modes of the proposed absorber in Figure 13. The 
power absorbed during arrest and the power expended during retraction 
are presented for a range of possible operating speeds. 

The maximum dimensionless power is ae shown as a function of 
Reynolds number Ry a Figure 14. This graph, although based on only 
three computed cee does correctly portray, approximately, how 


scaling affects power absorption. It would be a simple matter to 


calculate additional points on this curve from program NORMAL. 
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Figure 12 Non-dimensional Power vs. Non-dimensional 
Rotational Speed 
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As mentioned earlier, the actual existence of four stable modes is 
a surprising result. The original expectation, based largely on 
intuition rather than analysis, was that only one mode would be possible 
for each direction of shaft rotation, But the absorber is an unusual 
machine which tends to confound intuitions based on experience with more 
conventional devices. Consider for example the fact that losses are 
necessary for the machine to absorb the mechanical energy input. The 
smaller the loss coefficients, the higher the final equilibrium flow 
rate will be, and the greater the energy absorption. Thus we are led 
to the seeming paradox that the lower the loss coefficients the higher 
the energy lost! 

The power curves of Figures 12 and 13 are interesting and 
significant in that all curves have nearly constant and nearly equal 


slopes. The power can be expressed in the form 


where the slope of the curve is roughly 1\=3 for all cases, 

The horsepower curves of Figure 13 indicate that Mode 1, the 
design point mode, absorbs the greatest amount of energy. Mode 3, 
however, is a very close second. Refinement of the analysis should 
lower the predicted absorption of’Mode 3 somewhat. Nevertheless, the 
favorable comparison of this mode to the design mode is interesting. 
Figure 13 also shows that the power needed to retract the rotor is 
reasonably small for Modes 2 and 4, especialiy in comparison to the 
power absorbed during arrest. 

The eee of which mode the absorber will actually operate in 


under various conditions of use is interesting, but by no means 


completely clear at this time. It is thought possibly that if the 
absorber be started from rest with shaft rotation in the normal 
direction, Mode 1 swemmaroenut If retraction follows immediately, 

the retraction mode will probably be Mode 2, This is based on the 

fact that during arrest the flow will build up a large momentum which 
will tend to persist in the same sense during the subsequent retraction. 
If retraction be started after the fluid motion ceased, however, the 
resulting mode of operation remains uncertain. 

Some method of sensing and/or controlling flow direction would be 
useful to insure that Mode T was occurring on all arrests. Mode 3 
operation might not have enough absorption capacity to handle large 
aircraft. 

The general question of which mode occurs is an area for further 
study. Testing of the actual absorber with the idea of definitely 
determining the operating modes under various conditions would be 


extremely useful. 
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Ve. CONCLUSIONS AND RECOMMENDATIONS 


The development of the problem and the exact solution of the 
two sets of equations describing normal and reverse flow conditions 
clearly show that four stable operating modes exist. The analysis 
adequately predicts the absorption capacity of the proposed device 
for Mode 1. However, the analysis is not as accurate for the 
reverse flow or retraction modes. 

It is recommended that the retraction problem be further studied 
to refine loss coefficients and to more accurately predict the power 
necessary to retract the absorber. It is also recommended that 
further study be made into the question of predicting which fone 
occurs under various operating conditions. This study would 
encompass the unsteady phase of the problem, Actual testing of the 
absorber should also be conducted with the intention of determining 
which mode will occur. A method of men control or sensing would 
also be useful to guarantee proper performance in that Mode 3 
absorption capacity may be insufficient in certain cases. 

This turbo-type of energy absorber appears to be completely 
adequate to meet the demand of carrier operations now and in the 
forseeable future. The potential for future development and 


enlargement is ample. 
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APPENDIX A 


DETAILED ANALYSIS OF THE ENERGY. ABSORBER WITH THE 
FLUID CIRCULATING IN THE NORMAL DIRECTION, 
MODES 1 AND 2. 


A one dimensional mean streamline analysis after Vavra (Ref. 2) 
will be developed with the following assumptions: 


B 


constant 


2 


2 


Uniform flow at stations 1 to 4 except for blockage factors. 


constant 


Blockage factors are independent of the flow. 
The fluid is incompressible. 


Steady state conditions exist. 


1.0 Torgue and Velocity Relations 


The torque T, exerted on the fluid by the rotor is positive 


R 


when acting in the direction of rotation. 


 ceimoles Wie Vy] i (1) 
where 1 is the mass flow rate. 


m=p A ve 


p(2m R, bp) Yu, sas (2) 


p(27 R, bp) ‘ut, es 


The quantity 2n R b, is the flow area, i and Ve are the meridional 
velocities at the mean streamline, and Kay and Kno are blockage factors 


which are assumed to be independent of the flow. 
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By introducing a flow coefficient 


U, wR. 


all necessary velocity and thrust equations can be written in terms of 
known geometric parameters and the flow coefficient. An iteration 
technique is later employed in program NORMAL to calculate %. The 
following velocity relations are developed now in order to simplify 
the energy relationships which follow. 

The sign convention used is that all peripheral velocity com- 
ae are positive when in the direction of rotation, From Figures 


4 or 5 simple geometric relations are obtained which become, using 


Sunes. 
Wo = Vyp/eos B, = OR, é/cos B, (4) 
Wyo = Vy Tan B, = WR, & Tan B, (5) 
Vio = Wyo + OR, = OR, [1 + $ Tan 6] (6) 
Also 
Vo = Vyp Tan & = OR, % Tan a, We 


An expression for v, in terms of the known relative exit angle can be 


obtained by equating the previous two equations, or 


Tan @, = Tan 8, + 1/6 (8) 


similarly at the rotor inlet 


Vi aot (10) 


The mass flow equation and the definition of the flow 


coefficient yield 


R 
2 Keo 
V = 0 veo —_ C1) 
M1 2 R, Key 
and 
R 
2 Kp 
ut 2 R, Kay 1 1 


The mass flow equations must be extended to include flow 


through the stator. 


M27 Ry b, P Vigo Kao 277 Ay bs p Vung Ke 


(13) 


OTT Ry 1p 0 Vw Kea 


From the above relations, the meridional velocity equations are developed. 


“in . ea “pp 


ie oe (14) 
Yo "at “pe : 

Ty. Set Spe Ee 

Wp | cee | (15) 
Mg 0 fo Bae 2k | 3 
Vo Fa, “oy We ; (16) 


1.1 Frictional Losses in Upper Passage 


Let Mey be the frictional moment from forces along the walls 
of the upper passage from station 2 to 3, acting on the flow opposite 


to the direction of rotation. 


mR, Vi, = mR, Vip - Me (17) 
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Setting Mp, =A, DR V (18) 


and Tae = (1 - a) (19) 
eq. 17 becomes 
Ry 
Vg © a By Vue (20) 
From the mass flow equation 
Yu bon, = - - (21) 
3 eis 


jee  Bevetronal) losses) an lower Passage 


With M,, being the frictional moment caused by forces along 


£2 
the lower passage walls from station 4 to 1 acting on the flow opposite 
to the direction of rotation, the equations for the lower passage 


are developed in a manner similar to the previous development. 


Mog = hy m Ry vA (22) 
Toy = 1 Yg (23) 
R 
T aoapdiuca’: = hae (24) 
ut u4 Ry 
Hy TS b 
2 2 R 
vy, =? oR > oe = (25) 
MA 28, Ky by 


From Figure 4 it can be seen that 


vA = Vu Tan ah (26) 


Using the previous three equations gives 


R b 
RO as Goes a 


a Me oR Tan o (27) 


1 Ke = 4 


Di 


From Figures 4 or 5, one can relate the absolute inlet angle 


to absolute velocity components as 


Tan @, = =e (28) 


Using the relations developed for the inlet velocity eqs. 11 and 27, 


the absolute inlet angle is 


Je) 
R Sp 
Tan a, = — -=—— Tan @ (29) 
1° Tye B. a 4 
Similarly the relative inlet angle is 
V - wR 
Tan B, = os 1 
M1 
oR, 
Substituting for Vu yields 
b ie As 2 nies 
R 1 1 1 
Tan 8, = — — Tana, -(— — 1/6 (30) 
1 = "Me 3, et il qe) ae 


The relative and absolute velocities can now be obtained from 


their respective components by employing the following general equations: 


Z 2 2 
ia as (a) 
Z 2 zc ; 
Weevy +W, (32) 
W = V,/cos 8 : (33) 


De 


1.3 Torque Relations 


The torque as defined in eq. 1 can now be expressed in 
terms of known geometric parameters and the flow coefficient. Using 


ease 5(; 65 and 12 anveq, 7 yaelde 


B Ke 
pee a Ee pie “e404 {Tan Bo-TMhys oe Siete as (34) 


The non-dimensional torque coefficient +t, is now introduced as 


_ _ i ae 
Re > 5 b Des 
an(R,) = lo) Ko 
S 
fo) 
2 Keo R 
6+ 6° 4Tan B —— — Tana? (35) 
LTP Hosts Than ey, wha J 
The torque T which the stator exerts on the fluid is positive 
Ss 


when acting in the direction opposite to rotation 


4 mir gE) (36) 


Substituting the values for the whirl components gives 


Dd : ; me fee 
Sys 2 2 R Kno 
_ = en. z) p w ie $(1+9TanB,, ) - 0 \p_) ae Tan a 


Again introducing a non-dimensional torque coefficient 


a 
-_ : 
s Bee 
ae 2 
277 R a 
2 \z,) Kae Ce 
dD K, 
2 R 2 
= Thy 6 4 $ | Tym 22 Bo -— FF Tan o | 


pal 7 oy eit ¢ 


Met) Oy Se Stotveante ci<tern~s, wrvcd ‘to 


- a 


The general torque equation is 


The torque exerted on the rotor by the fluid minus frictional moment 
losses is equal to torque exerted on the fluid by the stator. 
The dimensionless torque coefficient introduced earlier is 


again given as 


Me 
= (38) 


2.0 Knergy Relations 


The torque and velocity equations developed in the preceeding 
section remain in terms of the unknown flow coefficient ®. In this 
section an algebraic energy equation will be formulated and solved 
to yield ¢. 

For the sake of this development, a total head term H will be 


defined which differs from the standard definition of total enthalpy 


En which is 
v- . 
Hy = b+ > (39) 
Ors 
2 BP 
Hy = ut +Deuts (40) 


where h is the static enthalpy 
p is the static pressure 
u is the internal energy 


and P is the total pressure. 
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Since the fluid is treated as incompressible, one can uncouple 
the thermal effects from the mechanical effects, and define a 'total 
head' H which does not include thermal effects. This new term is 


denoted by 


1 2 
eee vagiin 
0 Bz 2 (41) 


H=— 
p 


Obviously the enthalpies are related by 


(H, - u) =H (42) 


This new definition of total head H is useful in that constant H and 
constant Pa meiosis coincident. Throughout the remainder of the 
analysis, useful total head H will be referred to as simply head. Note 
- that this total head has the units of enthalpy (per unit mass). 

The peoneds undergone in the absorber is now represented in 
Figure 8. Obviously the process does not return to the original total 
enthalpy Hn4? sahara does return to the original total pressure 
and head. 

From Figure 8 the energy rise across the rotor is shown to be 


equal to the energy losses across the stator and through the passages, 


Ox. 


AH, - AH, - AHy - AH, = 0. (43) 


Therefore, 
(a, = H,) ny (H,, am H;) ee (a, 7% Hy) a (H, ay H,) = 0 (44) 


where all the above terms are written as positive quantities. 


DB 


m~ § 


alk ia rs: 


enPeelom lsemat! gheient Pon wood dot fot : 


The four terms of eq. 44 will be derived in terms of the flow 
coefficient and the known geometric parameters as in Section 1.0. Eq. 
hh is then straightforwardly solved to yield ¢. Program NORMAL is 
employed to accomplish this task. Since various loss parameters 
depend on $, an iteration technique must be used. 

2.1 Conditions in Rotor 
The incidence angle of the rotor blading i, is defined as 


R 


i= 6 fe  . (45) 


or merely the angle between the tangent to the mean camber line at the 
leading edge and the relative inlet velocity. The incidence angle i, 
is seen to become more negative in the direction corresponding to a 
greater lift coefficient for the blade. 

The relative velocity W, can be resolved into two components 
_respectively tangential and normal to the mean camber line. The energy 
associated with the tangential component W, is assumed to be largely 
recoverable. The energy pe cocueteed with the normal component is 


assumed to be totally lost. 


The useful relative velocity at the rotor inlet is then 


Wy = Wy ee 1p (46) 


and the velocity component lost due to incidence is similarly 


We 22 Woo sin 3 (47) 


The basic energy equation of turbo-machinery is 


ee oe ee eal (48) 


56 


which upon expansion becomes 


V, ry ai V. Eos a é 
z Ute eeimens We tie ee (49) 
Ho edie 2 2 
Using static enthalpies, regrouping, and cancelling terms in eq. 48 
gives 
We" - i Ww,” - me 
+ ——————— = _—_—_——— 
h, 5 h, + 5 (50) 


Substituting for the static enthalpy and rearranging yields 
2 2 2 2 


22 Rs ec ee aise Sa a (51) 
ij» 2)! es 2 2 al 


Upon resolving W,°/2 into its components and regrouping further, one 


obtains 
2 t 
p, wW Ato Wud we Ue woe 
(=++ -(=+5 + 2+) - ( a es | (52) 
5 ee 0 B 2 ees ae ? 


This leads naturally to the definition of the "equivalent relative 


total pressure" as 


2] 
gag 9) W. 
E2 ( 2 2] 
— = — + — 
; cB a ) (53) 
and 
12 Bei 2 
* EL = te ms + Be (54) 
p p 2 2 ? 


and to the definition of the rotor blading losses as 


7 2 
= 7 i 5 a c4 7 Ga ae, (55) 


It also proves useful to express this loss in terms of a dimensionless 


coefficient defined as 


- P 
fu 7 ase 
ies (CB ) (56) 
eae 


Recalling that 


ie = oo ae (42) 


and using the basic energy equation, eq. 48, one obtains the energy 


rise across the rotor in the form 
_ soe en 2 } 
(H, - HL) =U, Vo > Uy Vy - {sin ® aR oo oR ne (57) 


This equation can be conveniently non-dimensionalized by 
dividing by Uae The velocity expressions developed in Section 1.0 can 
not be utilized to simplify eq. 57, and the dimensionless energy rise 


across the rotor (IIA H). can be written as 


H, - H (cae it 


) eyak 
lun > vas Sel =) 


+ § {tan 3, - ea > Thy Ten o,) cos“ i = (58) 


nse 2 Boro 2 
PEI) +S Bw) 
+ *R 

2cos*B,, 
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where the dimensionless head coefficient is defined as 


Thi, pee (59) 


The loss coefficient Yp of eqs. 56 and 58 is a function of the 


deflection AS, where 


Ra ee 


and the aspect ration ps defined by 


n= (60) 


Tavs will: be further diseussed_in Section 3.0. 
2.2 Performance of Rotor Blading 
The rotor is designed for the maximum useful enthalpy rise 
across the blades. If this were a frictionless process, the useful 
energy rise would naturally be greater than in an actual process 
involving friction. Therefore, the rotor efficiency is defined as 
H, - H 


oe eae 


SS (61) 
on: 


“where the total heads Hy > H,» and Hom 


Figure 8. Of course the pressure rise through the rotor is related 


are shown on the H-s diagram, 


very simply to the corresponding energy rise according to the 


relation 
P 


i] = (62) 
p 


BY, 


The numerator of eq. 61 is obtained from either eqs. 57 or 
58. The denominator of eq. 61 is similarly obtained, but for the 
theoretical case both the incidence and blading losses are non-existent. 


It follows from eq. 58 that 


Il 
aN 


At. 


be Kae ee 
1+ 6 {tan Bo - Thee De Ky, Tan ay 


is the theoretical dimensionless energy rise across the rotor. The 


efficiency can obviously be written as 


(I AW R 5 
Tp = (ai, (64) 
fae 
The general definition of power ~ is 
(p=amdH (65) 
It is convenient to introduce a dimensionless power coefficient as 
a 
Ip (6 
PYviax Ro 
Recalling that 
V 
M2 
oo % (3) 
m= 27 p R, be Vip Kyo (2) 
and 
AH 
a0) R, 
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the dimensionless power and head coefficients can now be related as 


p= fon(G8) (67) 
The power is related to the torque by the simple expression 
p= The ( ‘ ) (68) 
MAX 


2.3 Energy Losses in the Upper Passage 
The energy loss in the upper passage, assumed as only a total 


pressure loss due to friction between the rotor and stator, is expressed 


as 
i P 
Te aug ‘ 
AH = Jol. él = —— i 6 
a D 3 5 : (69) 
or 
Py i P2 V5 
= (— + —) - (— + 
on iS 2 ) le 2 ) 7 (79) 
In a frictionless process 
2 
p V. 
— = ee > 
oe ray. ae 2 A) 
and from eq. 17 with the frictional moment Ma, equal to zero 
Ro ce 
Wa! 7"ue Rie (72) 
Spans re 


V 
a u3 
28 Pena 5 ia (73) 
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Similarly for an actual process 


2 2 
Pp ie V. 
rape 7) 


) 


the energy loss can now be written using the two previous equations as 


Pp = 12) ie 
Det, = 


2 ~ H, ; cam (75) 


a 
A convenient dimensionless pressure loss coefficient can now 


be denoted as 


Pp <1 1¢) 
C pee ee) . . (76) 


Gees 


70 
where oF will be shown to depend primarily on the radius ratio = of 


the upper passage in Figure 2. This is similar to energy loss 
‘coefficients developed for bends in smooth pipes. 


Eq. 75 is non-dimensionalized, and using eq. 76 one obtains 


° ° 
Cc UV Vv 
a ie dL gle Tae LUSe 
(1 AH) = > Esta G 1) ogee (77) 
— TW ee 


Using previously developed expressions for Vy one obtains an equation 


similar in form to eq. 58 or 


em 
ae 5 
a et 
ae 
2 
ie” Be 2 
Ss EN 2 e : (78) 
ASS ) tes (@) (Thy) Tan Bf 
nee fa Fo Pp Kyo2  (1-The) (ey enea 
—_—- — a, = an is 
ems EeR Se ae Rj as 


2.4 Performance of Stator Blading 


From Figure 4 the stator incidence angle i, is 


i = jor - @ B 
+5 3 3 (79) 
1 
and the useful inlet velocity We is 
! 
VW coske (80) 


The total head ahead of stator blades is given by 


2 
Ps Ve Pa 
J~ p < p ey) 
while the useful head is 
1,2 ! 
Pat Oe bee 
' 3 5) ue) 
H a 82 
D p 2 p (ee, 


The above three equations combine to yield 


(v,')° = (V5) 
3 = aE + 5 


yee (83) 
5 ZA 
HZ - ion (sin i.) 


3 


Similar to rotor blading losses, stator blading losses will be defined 


as 
! ! 
Peed Sk! ahs Mee (84) 
s Do 2 2 
yg | Male) 
or 
: 2 2 2 
Vv Vv Vv 
- qt! pars ine eA (85) 
Hy =H -Y,3- =H (sini e=ge itty SS 


Subtracting H, 


as - H, = H 


The same non-dimensionalizing technique used in Sections 2.1, 


is again employed and the energy loss in the lower passage is now 


available in the 


from both sides yields 


2 
Vv Vv 
ae ine eee 2s 
3 H, (sin 1) Daa Y, D 


general form as 


The above equation is actually 


He - H Ee - a 


ay ROOD har manele 
W Ry Ww Ry 


This can be further reduced to give 


which gives the energy losses in the lower passage in terms of the 
energy changes throughout the remainder of the system. 


step is to find the energy loss in the lower passage and solve eq. 


44, 


es eee 
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The next 


(86) 


anid) 255 


(88) 


Again the loss coefficient xe is a function of deflection and 


the aspect ration, where 


and 


2.5 Losses in Lower Passage 


The development of the energy loss in the lower passage 
closely follows the technique of Section 2.3, 


The energy loss in the lower passage is 


2 
p Vv 5 1 as 
p 2 


p 


for a frictionless process 


R 
ee: 


(Vaan > Sud R, 


or relating theoretical to actual velocities yields 


Therefore, 
2 
ce Ge ar en * wit ran 
Any ee p 2 23 
while 
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(89) 


(90) 


(91) 


(92) 


(93) 


(94) 


(95) 


the pressure loss coefficient is 


2 
ALE esta oo 
p Sa 
The expressions for Hy and A are now substituted into eq. 1, 
along with eq. 6 to yield 
é 2 
V. AB. 
e M1 2 = 2 
ST egres Wac ee C1 = Wy) \z,) Vus (97) 


This is easily put into the dimensionless form used previously or 


b Z 
(au), « 5 (g2) Te GE) + (1 = Myy?)(E BE nem oy) ] (98) 


Eq. 97 can now be substituted into eq. 88 and % can be solved 
for straightforwardly. Iteration is necessary, however, since all loss 
coefficients depend on 4, 

4.0 Loss Coefficient in Rotor 
Vavra (3) gives a relation for the losses in axial turbine 


bladings as 


2.28 4.97 
a ks baa cae ee aT lat oc (99) 


where AB is the flow deflection angle in degrees. The quantity ¥ is 
the velocity coefficient defined as the ratio of the actual and 


theoretical discharge velocities, which for the rotor is 


aan: (4100) 
Veo ees 100 
R Wom 
and 
AB, = 8, - B, (101) 
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The theoretical case occurs for YR equal zero which in 


accordance to eq. 56 gives P_, =P 


a ni ; : 
E1 ROTh’ Substituting eq. 100 into 


eq. 53 yields 


1 
Se Sap SP oes (102) 


The actual rotor blading losses using eq. 102 can now be written as 


Pn W : Dey Wes 
o if a2) = oa i) 1 
Yr = lames a! (103) 
2 ¥ 
| We /2 


In accordance with Horlock (Ref. 4) the loss coefficients depend 


on Reynolds numbers Rup and aspect ratio A, of eq. 60. Hence let Y. 


R R 
* 
>of eq. 103 now denoted as YR be the loss coefficient across the rotor 
5 
blading for Ap =~ -and Rip = 10°. The latter term is defined as 


ee | (104) 


where Dir is the hydraulic diameter of the discharge area of the flow 
channel between adjacent blades. IfL6 is the rotor blade spacings at. 


the radius Ros and ap is the opening at the discharge as shown in 


Figure 4, an approximation of the discharge opening is given by 


ap ~ A cos B,, j (105) 


then 


py ae Caen a oe (106) 
hR wetted perimeter 2 (a, + Dp) 
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Substituting for ap and rearranging yields 


2b, 


es = ae (107) 


1 eee cos 8 
2 
The blade spacing 4, can be found with respect to the radial 
extension Ry - Ry by using the blade loading criterion of Zweifel (Ref. ES 


For axial bladings this approximation is 


LO Pee a) 1 Miter = A oeaie (108) 
Ro % Ry (Tan Bo - Tan Bp) cos’B, 


Substituting this equation into the expression for the hydraulic 


diameter yields 


o( =) 
D \F 
a o Cae ae — Tan8,..)coss (109) 
ee BRIE? Geers Papaceee a 


0.425(1-R, /R,) 


Placing the above equation into the expression for the Reynolds number 


and recalling that 


oo (4) 


one obtains 


o(—*) 8 w ag 
0} | 


2 
(b,,/R,) (Lan8, - TanB,.) cos 2] . 


0.425 (1 - R,/R, 


Rup = 


(110) 


| cos8, + 


aa 


The loss coefficient is dependent on the aspect ratio and 


Reynolds number through the following relations after Horlock (Ref. Le 


Y, =—5-1 (111) 
2 Te 

t 1 2075 
Se ie CON Re ea (112) 


and 
AS, 


2 Tees 
ae G_) (113) 


! 
where the coefficient YR holds for the actual aspect ration only at a 


Reynolds number equal to 10°. By combining the above equations the 


generalized loss coefficient is simply 


- D225 
JA ena ' 207 (102 (aan 
XG = Te Cees | 3 3 - 1 <a Te ii 
ot (95 48) cn 


These relations allow the actual loss coefficient of the energy absorber 
to be determined under varied conditions. 
Obviously the Reynolds number is not totally specified by 
geometric parameters, but depends on the flow coefficient %. In 
1 


np? ora 


Reynolds number of 10°, and subsequently correct the loss coefficient 


iterating for ®, the process must be assumed to begin at Y 


as improved values of 3 and Rup are calculated. This procedure must 
also be employed for the velocity coefficient ¥ which is dependent on 
$6 through a dependence on the actual inlet angle Bae The iteration 
must be continued until convergence is achieved for all values 


dependent on $, 
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Sette amitifales poly) 


3-1 Loss Coefficient in Stator 
This section follows closely the development in Section 43.0, 


Eq. 99 is again used with 


Yoere (115) 
4Th 
and 
AB. = el =a 6 
Ba = 0s men (116) 
For a frictionless process in the stator, the total pressure 
ae 
P : i ; 
(Pagan at stator exit would equal total useful inlet pressure Png ; 
or 


i . 
Coon _ Poa _ Pa , Ma) on (117) 
(ine Oe 0 2 


2 
12 p V 
eels Se ee) gas (118) 
p p y D) 
Ss 
For an actual flow 
2 : 
ie) p V 
es coe en (119) 
a a ns ee 
which yields through eq. 84 
1 2 
= -2 atl aad (120) 
Ss 


This was as expected from Section 3,0. The loss coefficient can now be 


written after eqs. 112, 113, and 114 as 


nee £07 | ie 
ae ee oe can 
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402\" 29 


reas, Ge) (122) 
and finally 
ig HOB le eee : 
ies ye (915 + mM ey ) 2) 


where the aspect ratio is defined by eq. 90 as 


Die Ve 
hs 4 
Ry a = (124) 
Ss v ; 
where for a. nee oe the hydraulic diameter is 
Ale 
D a be eee (492 
hs b [a (125) 
4 sis 
cos @ 


A 


The blade spacing according to Zweifel (Ref. 5) is 


Ssoanaen 0.425 


eon. j 2 (126) 
3 fee | ar oie Tan a, |)cos 2 
and from eqs. 25 and 26 
RR, P K, 
2 R 2 1 
Ve«twuRe- Soe o_o (27) 
4 2 “i be. Koa cos a 


The Reynolds number can now be expressed using the above four equations 


as 


fj 


o 


-D 
aaa ta 


So AG A ee 
, es 
cos G, + 5 = il Tan a, - Tan a,, |) cos ay (128) 
R 
4 
0.425 (R,/25) (1 - x) 


A series of iterations is again needed to insure proper 


convergence of all parameters, 


4.0 Pressure Loss Coefficients in Upper and Lower Passages 


In accordance with eq. 76 the pressure losses in the upper 


passage are expressed by 


epee a2 = OC I= (129) 


and similarly for the lower passage by eq. 96 as 


2 
(o,) =p V. 
(iedat 1 M1 


With Apb being the pressure drop in the bend, the loss 
coefficients of these bends defined as 
~ 2 y2 (131) 


are shown in Figure 15 as functions of the Reynolds number Rys where 


similar to eqs. 104 and 124 


Ry = 2 (432) 


Vavra (Ref. 2) has developed the set of curves shown in Figure 15. 


These curves are based on the data of Sprenger (Ref. 6). Interpolation 


fe. 


was to obtain values for ratios of r/o different than those found 


in Sprenger, 


The hydraulic diameter Dd is again 


4 flow area 
Dae (133) 
wetted perimeter 


For a ring channel with flow area of 27 R b the wetted perimeter is 
2(2nR), ox 


on R)b 
Bie = ca 7 = 2b (134) 


Hence at the rotor inlet there is 


D = 2b, (135) 
and at the stator inlet 
‘ DE = 2b. (136) 


These above two values are introduced into eq. 132, with V 


equal to the corresponding value of Vu 


The result is 


that is, Vu before the rotor — 


and Vu before the stator. 


V. 2b 
M 
R= (137) 
V. 2b 
M1 R 
ee (138) 


The Reynolds number is then used to obtain the loss coefficients §€ of 


Figure 15. 


The loss coefficients € of Figure 15 have been obtained f 


measured pressure drops by subratcting the pressure drops that 


12) 


“s 


occur in a straight duct at the corresponding Reynolds number and 
surface roughness. 

There are several conclusions which can be drawn from Figure 15 
and from additional data in Sprenger (Ref. 6). In general, the loss 
coefficients for bends made up of two concentric radii increase with 
decreasing Ti/b ratios. The curves€ versus Ry that exhibit a 
minimum indicate the flow separations occur at the inner radius, and 
should be avoided. Finally, Sprenger (Ref. 6) shows that two bends 
with Ti/b = 0.5 arranged in series to form a 180° bend have a loss 
coefficient of 1.66§ at Ry = 10° where € is the loss coefficient 
of one of the 90° bends. This data will be a first approximation in 
choosing the loss coefficient for the 180° bends in the energy 


absorber. 


and Lower Passage 


The losses in the upper and lower passages will be assumed 
due to two components, aoesaae pasasne and TeLetLOl. The losses due 
to duct bending were discussed in the previous section. An arbitrary 
assumption is to take these loss coefficients equal to 2.0 and 2.5 
for the upper and lower passage respectively. The values of € are 
given in Figure 15 for various value of *i/ be It is recommended that 
Ti/b be chosen between 0.5 and 1.0. Any bend with a value of ri/b 
larger than 1.0 will be assumed to have the same losses as a bend with 
Ti/d equal 1.0. 

To the above losses must be added the frictional losses in a 


straight duct of hydraulic diameter DO and average length L, where 


Tee oe (139) 


; ia 
a we é Ves fT eT | 7 alae 
1? ss) | a 


. 


and L; is the distance along the inner contour and Ly that along the 
outer contour in the meridional plane from the trailing edges of one 
blade row to the leading edges of the other row. (See Figure 2) 

In accordance with eqs. 132 and 136 this general frictional 


pressure loss is 


= ie (140) 


The friction factor }, is a function of both Reynolds number 


ai 


and surface roughness. The value of hp is obtained from the curve of 


Figure 16, also from Vavra (Ref. 2), for a relative surface roughness 


of 
K 
na i 
7 ae (141) 
ee 10 


where K. is the peak-to-valley roughness of the surface. Hence, for 


this case 


2b | + 
K. = TOF = 200 b u. = fin. (142) 


If it is assumed that the rms roughness be one-half Ks the surface 
finish must be better than 100 b micro-inches. 

The applicability of this data to the three dimensional flow 
channels of the energy absorbed is obviously questionable. The losses 
defined below must , at best, be considered as first estimates only. 


The loss coefficient c. of eq. 129 will be taken as 


L 
2 u 
C = 2.08 + 2b he (143) 
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"3 shal d 


OLOVS NOG Jid3 


(2 *yoy eanre)) Jy xoqoeg UOTLOTT TTEM Teuueyp peyeuTysH 9) emmsTg” 


NY YSIWNN SGIONASY 


8 Fe G «ve c ¢e Ol 8 2= Ss v . 3] 


a 


and Cy of jeg ua30 dis 


C= 203 eer (144) 


where € is obtained from Figure 15 and hp from Figure 16 for 


2p. V R bee 25 
s ‘M3 Bie BO aOR NE be 
SS oO oR, (145) 
Ryu v ( ] R 3 =) v 
and 
2b_ V R op 
R ‘M1 2 ‘po 
Ruy = = (g0R, = *B2 _ (146) 


pee 


In eq. 4 the length L, is given by 


(Coes (Cae) 
Li = aa WL (©) 


where (a), and age are the lengths of the inner and outer 
meridional contours, respectively, of the walls of the upper passage 


between stations 2 and 3. Similarly in eq. 133 


ee gt le) 
wera: ho 
Ly = 5 | (147) 
where (Ly); and (Ly), are the counterparts of (ene and (Lo for the 
lower passage between stations 4 and l. 
5.0 Losses of Angular Momentum in Upper and Lower Passages 
The frictional moments along the upper and lower passages 


were introduced by eqs. 18 and 22 respectively as 


hy a es | (148) 
and 
= ’ ho 
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No data is presently available for the losses of angular 
momentum in passages of the type needed between the rotor and stator 
of the present absorber. A very crude approximation of these moments 
is obtained, however, by assuming that the shear stresses se on the 
walls are merely those acting on a flat plate with a peripheral velocity 
Me at R being the velocity outside the boundary layer. This over- 
simplification ignores the effect of the meridional velocity components 
which may cause large increases in the frictional moments. 

For fully turbulent boundary layers along rough walls the 


shear stress is 


where Cp is assumed to be invarient with the Reynolds number of the 
fiow. On the element of wall area 2n Rd L of Figure 2 there acts a 


frictional moment 


dM G27 Rolin Vo eR 
af Oo 


which becomes when using eq. 3 


It is assumed that the passages have been designed such that the flow 
is not separated from the walls, and all stresses external to the 
boundary layer are an order of magnitude smaller than stresses along 
the walls, and therefore negligible. This results in the product (RV,,) 


being constant outside of the boundary layers and equal to R the 


2Vue? 
value at the rotor discharge, for the upper passage. Integrating 


eq. 151 from station 2 to station 3 yields 


te 


Mie weg" eg? Apew ot hefect 


eieris “oss A 


(3) 


Mee = sles Tipe Ee ees dl, (152) 
Ha = 7p es Re Wa Les aia). (153) 


where the inner wall is of length ae and the outer wall is of 


length Ci) Hence by eq. 148 


[y)y + Oy) 
— UN ee Lee oe = (15!) 


m 


Substitution of the mass flow equation and velocity relations into 


eq. 154 gives 


Oe TE Bo [a.), % (1), 


SET Sina it ee EET (155) 
6 (2 dp Ky 
from eq. 19 
r 
1+ 3 Tan B (Ch. Gin) | 
Coase ig et u’o 
Me Gye eee eS ee Lee (156) 
Thytu f 5 20s Kes 


For the lower passage a similar procedure and derivation leads 
to 
- 
+ 
ee Tan a, | (Ly); (Ly), | (sr) 
coe Pow 


oh s 


and 
Cp, Tan m, [ (4); af (Ly), | 
2] Ka b 


Thyp = ie (158) 
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In eqs. 155 through 158, the absolute value of the flow angles 
must be introduced since the loss is a function of sbsolute deflection. 
For fully turbulent boundary layers, as in the present case, 


the local skin friction coefficients c, depend on the surface roughness 


fe 


ration K/%s where K. is the peak to average value roughness and x is 


length along the plate. 


ec, = .OO4 for K/x 10 


1073 


Cp = .006 for K/x 


Since there wiil be flow irregularities and possible flow-separations, 
go ABS .016 was chosen 


since it gives best agreement with M-21 test data. 


it is advisable to use a larger value for c 
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APPENDIX B 


DETATLED ANALYSIS OF THE ENERGY ABSORBER WITH 
THE FLUID CIRCULATING IN THE REVERSE DIRECTION, MODES 3 AND 4 


A separate set of equations must be derived for the reverse flow 
modes. The reason, mentioned earlier, is that the relative flow exit 
angles which are assumed constant are now different. 

This development follows very closely to the previous analysis. 
Therefore, the intermediate steps of this derivation are omitted. 


1. Assumptions 


a. By = constant 
a = constant 


The remaining assumptions are identical to the normal flow case. 
2. Analysis 
The geometric relations employed are obtainable from Figures 6 


' 
and ‘7. A new flow coefficient ¢ is conveniently defined as 


eee ae (1) 


The torque exerted on the fluid by the rotor when acting in the 


ined 


direction of: rotation is 


meh Sehr athe Ha SS) (2) 


where the mass flow equations are identical to those for the normal 


case. 
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The frictional moments acting from stations 1 to 4 and stations 


3 to 2 are respectively 


FL 


Ss 
2 

I 
a 

Be 
Hs 

< 
oN) 


where 


Tk yu fee 
m Ro Vie = ii R3 oe - Mey 


The relative flow inlet angles are 


and 
(1 + 6° Tan Be) 


Tan. Q = -—O_ 
aa. enter 


yu Ps 
The torque is now reduced to 


b 
» °() 2 
Tp =27m1 0 R, Ry Ka Ww Tp 


where 


84 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


Proceeding to the rotor, one can write, 


! 


Wy = Wo cos ip 


ee ees = 
Pro y Pay 

YE = 

R ae 
exer 


and finally 


1 2 he 2 R =) 
pan? wR © . Ae 
w R al 
i: 
b 
SOR 5 
+ 6 {tan B. - (= — Tan @ ) cos aL } 
i KR3 Be "wu 3 R (9) 
. Diy —_ 2 r nN aa 7 in LY 
ee sah ay Bl 31 °R a *R 
8 { 5 a) i ) + Gay me) |+ =—} 
2 Keo a3 2cos By 
where 
wa is 
os 
R, - Ry 
for this case 
AH x 
WAVHI 2B (10) 
w R 


In the lower passage 


p = 18 
Sth 4 (21) 
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and the energy drop is 


2 

WS Eh ie wg) (1. 
a a 2 R 

wR 4 


8 {01 = my )Ggt) ae 04} 


é Dee 2. le n° ) A 
ee =) eG) ey 


Now fom) the stauore 


Vy = Vy cos i, 
ve =e pes A 
Poa 7 Pog 
Ne = = 
s ey 
2» 
ee 
A = 
R, -—R 
Suae A 
A (ose es (04 
Behar sr= 8p 


and the energy drop across the stator is 


68 Gy) GE eB) ome 


R 


= ; 2 us ae b_ a) 


(12) 


(13) 


In the upper passage 


p ria ©) 
Ge 2Th 2 


The loss coefficients 


CL = 2,0 5. —— 
Cy = 26S + 

Ay = Ce an @ 

hg = on Tan B 


are now given as 


The Reynolds numbers ot obtain & and h 


fy eee ae ae 


i Re 


The expressions for Y, Yoo 


analysis, but the Reynolds numbers are now 
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graphically are 


(14) 


(15) 


(16) 


a) 


(18) 


and Me remain the same as in the previous 


! 2 
w R, 


be <8 
a v 
1 
(b_/R | Tan Bi ete 16 | eos B 
Rout 1 a2 1 


0.425 = 1) 


Rup = 


cosB) + 


(19) 


ey 


2 | Tan OB - Tan a | cos‘, 
ase + 1 
Bay (83 
0.425 (e) x - 1.0) 


Rs = 


The equation which is finally solved to obtain 6 is 


A Hp a q + A H, + A Hy 
Fox 
i 1 (Hy - H,) - (H, = H) = (H, _ H)) = 0 (20) 


The equation is solved by the program. An iteration technique must be 


! 
se 


carried out wntil all parameters dependent on ® converge. 
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APPENDIX C 


PROGRAM NORMAL 


Program NORMAL is essentially a modification of the program 
developed by NAEKC. 

Program NORMAL solves the equations governing normal circulation, 
Modes 1 and 2, for a given design. The solution is based on an 
iteration technique which is repeated until the flow coefficient 
$ converges to within a specific accuracy. A starting value of $ 
is picked as are initial values of the various Reynolds numbers which 
are also ® dependent. The initial iteration step occurs with all 
losses based on standard initial Reynolds numbers. This subsequent 


value of $ is then used to improve the estimates of the Reynolds 


¢ 


numbers and consequential losses, The technique is carried out 
until ¢ converges to the proper value. 

If the design parameters and angular speed are input data, the 
program will compute the flow coefficient, torque, velocities, losses, 
energy changes over each component, and the dimensionless torque, 
power, speed, and energy. The flow chart of Table IV depicts this 
process, 

The retied of data input is clearly indicated in the program 
itself, — The program is felt to be self-explanatory in that comment 
cards have been extensively used. The Fortran symbols employed are 
given in the list of symbols. 

Subroutines have been used liberally. Each subroutine accomplishes 


one individual task such as determining a loss coefficient, or 
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enthalpy change across the stator, etc. The subroutines Fig. 15 and 
Fig. 16 are closest point approximations of the curves depicted in 
Figures 15 and 16 of the analysis. Thses points are given by Block 


data which is called by each subroutine. 
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Table IV SLOW CHART FOR COMPIITER PROGRAMS NORMAL AND REVERSE 


Data Input: Desien size - 


Assien initial values. of 
$ and Reynolds numbers 
to begin iteration, 


New ¢@ and 
Revnolds Nos 
now used. 


a 


Tan =C > D 


sau Cag 


2A ee Oe ee 
TEST, N= | 


| PASS TEST 


Flow Analysis: Calculates and Prints 


Enercy Differences 

Losses and Loss Coefficients 
Non-dimensionalizes all pertinent quantities 
Pressure Differences 

Torque 

Flow Coefficient 


{ 
Velocities | 
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APPENDIX D 


‘PROGRAM REVERSE 


Program REVERSE solves the equations for reverse circulation, 
Modes 3 and 4. The program is very similar to program NORMAL, and 
the flow chart of Table IV is the same. Further discussion would be 


mere repetition of Appendix C. 
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